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ABSTRACT 


Measurements of the reflected solar radiation and 
emitted infrared radiation were obtained over the entire 
globe from the Nimbus U meteorological satellite during 
the period from 16 May 1966 to 28 July 1966, From these 
measurements, the outgoing long-wave radiation flux* the 
albedo* and the radiation balance of the earth- atmosphere 
system were computed for five subperiods* each of a half 
month’s length* All results are presented here in maps* 

The global albedo was found to be between 29 and 31 
percent* which is less than earlier accepted values of 33 
percent and more. The globally emitted long -wave radia- 
tion flux results in an equivalent black-body temperature 
of the planet earth of about 256°K* which is 3—4 degrees 
higher than was found in earlier studies. These results 
indicate that in earlier studies of the radiation balance 
the cloud cover and possibly the cloud height were slightly 
overestimated. The global averages of the radiation bal- 
ance obtained were between 0*002 cal cm* 2 min* 1 of surplus 
and 0,007 cal cmT 2 min -1 deficit* 
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THE RADIATION BALANCE OF THE EARTH-ATMOSPHERE 
SYSTEM FROM RADIATION MEASUREMENTS 
OF THE NIMBUS II METEOROLOGICAL SATELLITE* 

by 

Ehrhard Raschke 1, 

Goddard Space Flight Center 


INTRODUCTION 

In studies of the heat -balance of the earth- atmosphere system* only two energy fluxes are con- 
sidered: radiative energy of absorbed solar electromagnetic radiation and of re-emitted long- wave 
(thermal) radiation. Other possible energy sources and sinks (such as heat flow from the earth's 
interior and crust and solar particle radiation* as well as energy transport by waves to space) are 
unimportant and can be neglected here. 

About 99 percent of the energy of solar electromagnetic radiation reaching the top of the at- 
mosphere is within the spectral range between 0.2 and 4.0 microns. Its shortwave component* for 
wavelengths <0*3 microns (which amounts to 1.2 percent of the total flux)* is almost completely ab- 
sorbed in atmospheric layers 50 km above the earth’s surface. Of the rest, about 60—70 percent 
are absorbed mainly at the ground and in tropospheric layers, while 30—40 percent are scattered 
and reflected from ground to space (Reference 1), 

Of the flux of outgoing long-wave radiation, only about 5 percent originate in layers above the 
tropopause; the main part is emitted from the ground and from tropospheric layers to space. Thus, 
any investigations of the radiative budget or radiation balance of the earth -atmosphere system will 
consider mainly the radiation balance at a level less than 30 to 50 km above the surface. Results on 
spatial and temporal changes will be highly correlated with properties of the system beneath this level. 

Our present knowledge of the radiation balance on the earth- atmosphere system was obtained 
from earlier investigations* whose source material was climatological data (References 1—7), 

These showed that at low latitudes the absorbed radiation from the sun exceeds the re- emitted 
thermal radiation lost to space, while over both polar regions more radiation is lost to space than 
is gained by absorption (References 8* 9* 10), The thermal gradient between these regions of energy 
surplus and deficit causes the circulation of our atmosphere. Since the earth- atmosphere mean 
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temperature did not change over periods of a few years, it was assumed to be at radiative equilib- 
rium at these times. 

The accuracy of the earlier investigations was limited by the sparse material for observation 
available at that time; this affected their usefulness in the study of the atmospheric circulation. 
Therefore, satellite experiments were designed to measure the reflected solar radiation and emitted 
infrared radiation continuously over the entire globe. Such experiments had already been flown 
aboard Explorer VU {References 11, 12, 13), and several TIROS satellites {References 14 through 
18), but these satellites were not in polar orbits; their measurements covered only parts of the 
globe between about 60°N and 60°S. Reliable results were only obtained for the outgoing long- wave 
radiation. Results obtained for the reflected component of incident solar radiation showed large 
uncertainties. These can be partly explained by postlaunch degradation of instrumental response 
{Reference 14} and by the assumption of complete isotropic reflection characteristics {Lumber t- 
surface) of the earth- atmosphere system. There is so far no report on results obtained from meas- 
urements of the Soviet COSMOS- satellites. 

The meteorological satellite Nimbus H (Reference 19) was the first satellite to measure the 
radiation balance over the entire globe. Launched on 15 May 1966, its instrumentation made meas- 
urements from 16 May to 28 July 1966 (a period centered on the northern summer solstice). From 
these measurements were computed the outgoing long- wave radiation flux, the reflected solar- 
radiation flux, and the radiation balance. The solar constant was assumed to 2.0 cal cm" 2 min" 1 . 

The anisotropic nature of the earth- atmosphere reflection characteristic was taken into considera- 
tion, with simple models derived empirically from published data of airplane and balloon measure- 
ments. The computational procedures are discussed extensively in the next three sections of this 
report Appendixes A, B, and C present some additional formulas used in the computer programs, 
and describe the determination of empirical models on the anisotropy of the reflection properties of 
the earth- atmosphere system. Appendix D summarizes the data flow in the computer program. 

Results on the albedo, the outgoing long -wave radiation* and the radiation balance of the earth - 
atmosphere system were obtained for five semimonthly periods covering the full period of available 
Nimbus H measurements. Appendix E shows maps of these quantities; they are briefly described 
in the main text under '’Geographic Distribution of Results,’’ Geographic distributions for selected 
periods have received a more detailed discussion elsewhere (References 20 and 21) and will be the 
subject of a forthcoming paper by Raschke and Bandeem 

The global and hemispherical budget is discussed further on, while the subsequent sections 
make comparisons of the new results with earlier investigations and with direct measurements of 
emitted thermal and reflected solar radiation. 


THE NIMBUS II MEASUREMENTS 

Nimbus n crossed the equator on a nearly polar, sunsynchronous, circular orbit, northbound 
near local noon and southbound near local midnight; its mean height above the earth's surface was 
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about 1140 km. Because of its orbital period of 108.6 minutes, the entire globe could be observed, 
day and night, within a 24-hour interval. The reflected solar radiation and the outgoing long-wave 
radiation were measured in the spectral ranges from 0.2 to 4*0 microns and from 5,0 to 30,0 microns 
with a medium -resolution radiometer, which is described in detail elsewhere (Reference 22). 

The field of view of the scanning mean- resolutioo-infrared- radiometer was about 2,5 degrees. It 
enabled a spatial resolution of 50 km (near the subsatellite point) to about 110 km (at a nadir angle 
of 40 degrees). The radiometer scans were performed perpendicular to the orbital plane. The 
measurements of outgoing infrared radiation were continuously controlled by an onboard calibrated 
system; thus they are accurate within ±2 percent. The consistency of measurements of reflected 
solar radiation was investigated, comparing the bidirectional reflectances (Equation 3, page 6) over 
various cloud- free portions of the Sahara and Saudi Arabia, That is, the ratio between the measured 
radiance and the vertically incident ir radiance of solar radiation was determined {see next section). 
These cloud- free portions were selected according to the equivalent black- body temperature of 
concurrent infrared-radiation measurements in the 16- to 11-micron range. The bidirectional re- 
flectance over these areas varied within limits of ±5 percent of a mean value, not indicating a tem- 
porary change of instrumental response. Thus the measurements of reflected solar radiation 
should be fairly accurate, provided that there was no deterioration of instrumental response im- 
mediately after launch. 


THE BALANCE EQUATION 

The net radiation flux, Q, or radiation balance of the earth-atmosphere system at a horizontal 
surface element outside the atmosphere of geographic longitude, K , and latitude, 4 > , is the sum of 
radiation fluxes crossing this surface: 


Q(\, 4>, d) = S(A., <t>, d) - R(\, <t>, d) - E(\, d>, d) , (1) 

where S and R are the fluxes of incoming and reflected solar radiation, respectively, and E is the 
emitted long-wave radiation flux. The letter d designates the day of the year for which the balance 
computations are made. {It stands for the sun's declination, the true distance from sun to earth, 
and the measurements at this particular day “according to the subject under discussion.) The dif- 
ference s-R determines that part of the sun's radiation which is absorbed in the earth-atmosphere 
system, while the ratio r/s is the albedo. 

The computation of the outgoing fluxes from measured radiances requires three major steps: 

1. Conversion of the measured filtered radiance N ( to unfiltered radiance N, 

2. Integration over all angles of measurement, and 

3. Averaging over a 24 -hour interval to obtain mutually comparable values. 

Since the earth- atmosphere system up to only 30 to 50 km above the surface is almost the only 
source of reflected and emitted outgoing radiation fluxes R and E, the second step can be simplified. 
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The radiation from the earth -atmosphere sys- 
tem is considered as coming entirely from the 
bottom of the atmosphere, i,e,, from the earth's 
surface (a sphere 6371 km in radius). Then, 
the integration of outgoing radiances n will be 
carried out above an observed surface element 
by 

J * 2 7? fir/2 

| N(# , sin # cos 9 d# , (2) 

o 1 o 

where $ and \p are the zenith and azimuthal 
angle of measurement respectively {Figure 1). 
Uncertainties caused by this and other simpli- 
fications will be diminished sufficiently by av- 
eraging over all single results falling within 
gridfields about 300 km X 300 km. This grid size has been used in mapping the results 
(see Appendix E), 

COMPUTATION OF THE FLUX OF 
OUTGOING LONG-WAVE RADIATION E 

The outgoing long-wave radiation E was computed as shown in Reference 23. The first step 
was to convert the measured filtered radiance to the unfiltered radiance (Le, ? the infrared radia- 
tion emitted in the spectral range from about 4 to 50 microns); Lieneseh* derived a relationship 
between radiances N f and N for about 60 atmospheric profiles. Figure 2 shows his results for a 
vertical upward path (plain dots) and for a slant path (circled dots) with a zenith angle of $ = 78,5°. 
These results show that the relation between N f and N, is nearly linear and independent of the zenith 
angle 6, since the spectral response of the Nimbus n Medium Resolution Infrared Radiometer (MRIR) 
{5,0 to 30.0 microns) covers a region where 80 percent of the total radiative energy is emitted to space. 

The further integration of n was accomplished assuming symmetrical radiant emittance at the 
viewed area. Thus, only the dependence on the zenith angle 9 needs to be known. It was obtained 
statistically by Lieneseh and Wark (Reference 24) from measurements of TIROS satellites. The 
formulas are presented in Appendix B. 

The computation of a 24-hour average, as expressed in the third of the above steps, requires 
several measurements over the same area at different local times. Diurnal variation of atmospheric 
conditions (mainly clouds and temperature) over some areas causes diurnal variation of the outgoing 


zenith 



Figure 1 — Geometrical configuration 
of earth and safe! life. 


*Lienescb, f - H. , pr i v,* ;c comEUiiaicadof], 1966 , 
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Figure 2— Relation between the filtered (5.0—30,0 microns) and the 
unf titered radiance of outgoing long-wave radiation. 


long- wave radiation. The Nimbus II measurements indicated that nearly all land areas were con- 
siderably "warmer 11 by day than at night* while over the oceans the difference between day and 
night measurements was in either direction. 

But between 40 °N and 40 °S, the satellite orbit permitted only one or two measurements by day* 
near local noon* and one or two at night near local midnight. More measurements were obtained over both 
polar regions, where the sufeorbital paths overlapped each other. These measurements, were as- 
sumed to be representative for each condition, day and night The 24- hour average* then* was ob- 
tained by weighting each flux value over the same gridfield with the length of day or night over this 
particular area. This was done with the semi-monthly averages of flux values for each gridfield. 
Over areas where either day or night data were missing, no values of the outgoing flux were de- 
termined, Thus there are gaps in the geographical distributions of the outgoing long- wave radiation 
flux. They are indicated in the maps by dashed xsolines and stippled areas (Appendix B). 

COMPUTATION OF THE INCIDENT [S] 

AND REFLECTED [R] FLUX OF SOLAR RADIATION 

Theoretical (References 25 through 29) and experimental (References 30 -32) studies of the 
reflection characteristics of the earth -atmosphere system showed a strong dependence on 8 f 
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and £ (Figure 1), This was confirmed by statistical analyses of satellite measurements of reflected 
solar radiation Levine* and References 33 through 35, Viezee et al. (Reference 36) and Rabbe 
(Reference 37) demonstrated the angular dependence of reflected solar radiation, with examples of 
satellite measurements. 

Thus, the bidirectional reflectance p ( ' of an observed area K, 4> (as obtained from N f , the meas- 
ured radiance of reflected solar radiation) must be considered dependent on 6, </>, and £' 

N f (\, <t > ) 

P s ' 'P, £' \ *•, 4>) ~ 77Tc , steradian" 1 , (3) 

COS ty ' Of 

where S f ' is the incident solar irradiance at the moment of measurement and within the filter range 
of the instrument. It is obtained by integrating over all wavelengths the product of the spectral 
response of the instrument (Reference 22) with Johnson's data (Reference 38) on the extraterres- 
trial spectral solar irradiance. The latter are based on the solar constant s 0 = 2.0 cal cm" 2 min" 1 . 
Then, the solar irradiance on the day of measurement was determined by 

S' = s 0 (f 2 /f 2 ) (4) 

where l 1 and t are the mean and the true distance from sun to earth, respectively. 

The response of the instrument covers a spectral range, within which about 99 percent of the 
energy of solar electromagnetic radiation enters the atmosphere. Therefore, it is assumed in all 
further steps that p f * is a mean value for the entire solar spectrum ( p { # - p * ). As for outgoing 
flux of emitted long-wave radiation, the computation of the reflected flux of solar radiation requires 
measurements under various values of <9, 0, and V , However, over most areas of the globe only 
one measurement was obtained in 24 hours. Thus, it was required to determine the outgoing flux 
from a single radiance measurement. The steps in this process are discussed below. 

We compute the directional reflectance r {the ratio between the solar radiation reflected into 
the upper hemisphere and the incoming solar radiation} from the bidirectional reflectance (Equa- 
tion 3), using empirically determined relations. These relations were derived from airplane and 
balloon measurements of reflected solar radiation with different values of 6 and ^ (References 30 - 
32). Figure 3 shows one of three diagrams in which isolines X = r/p (where p - v * p') are drawn 
vs 8 and This figure was used to convert all p f from measurements at solar zenith angles £', 
such that 60° < V S 80°, to r. Two other diagrams were constructed forO < C* - 35° and 35° < V - 60° 
(see Appendix C). The patterns of the isolines X a r/p in Figure 3 shows that at low sun an ob- 
served area appears much brighter near the horizon than it would be at the subobserved point. The 
integral of r/p over the hemisphere is The outgoing flux R f of solar radiation reflected from a 


*Levine f j. S* t “The Planetary Albedo Based on Satellite Measurements Taking into Account the Anisotropic Nature of the Reflected 
and Backscattered Solar Radiation,” Master’s Thesis, New York Univ., Graduate School of Arts and Science (private communication). 
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surface element to space between sunrise t r 
and sunset t is then obtained from 


*'<*-.*> = r -4 T - I* * r(e(t);\,4) 

s r J 


S' cos £(t) dt , (5) 


180 ° 

0 = 90 * 



0= 30” 

X ( 0. CONST. ) i 

1 .55 

r ( C = CONST. > \ 

p {$,*,<-■ CONST.) ; 


0=30° 




0 - 60* 


e -- 90° 

0-J- 
- 0 


J0* 


0 = 0° 


>= 135 ° 


I = 90° 


0 - 45 ° 


where t s - 1 r is the daylight period in hours of 
an observed area (k, 4 )» Here a general rela- 
tion r(£) was required* which describes the 
change of the directional reflectance r with the 
sun’s zenith angle; this relation was derived 
from the above measurements and from re- 
sults of statistical analyses of satellites meas- 
urements (References 33 and Levine*)* Fig- 
ure 4 shows a normalized plot r (£)A(E - 0) 
vs £ of this relation and compares it with the 

statistical results of other authors* The directional reflectance of a surface element was found to 
increase with increase of the sun’s zenith angle* The 24 -hour average of the outgoing flux of re- 
flected solar radiation is: 


Figure 3— Dependence of the ratio X = r/p on the angles 
0 and of measurement, at very low sun (60° < £ S 80°)* 


R = (t,-t r ) • R'/24 . (6) 

Appendix C explains the determination of the above-mentioned models in greater detail* 

Using the relation of r(£)A(£ = 0) to £ in Figure 4 to compute the flux of solar radiation re- 
flected from an area to space between sunrise and sunset implies this assumption: that the physical 
properties (cloudiness and surface state) of an observed area do not change throughout the day* 

This assumption introduces a major error into the results over areas with a diurnal cloud variance 
that is stronger than the day-to-day variance due to transient disturbances* This error is most 
likely to occur over low -latitudes* where only one measurement was obtained per area per day* 

The incoming fluxes of solar radiation can also be computed from Equation 5, setting r = 1. 

In integrating the equation, refraction in the atmosphere was taken into account with a simple model 
derived from Table IX in Reference 39* 


•Levine, J* S*, "The Planetary Albedo Based on Satellite Measurements Taking into Account the Anisotropic Nature of the Reflected 
and Backscattered Solar Radiation,* Master*® Thesis, New York Univ*, Graduate School of Arts and Science (private communication)* 
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Figure 4— The relative change of the directional reflectance (r) 
of the earth -atmosphere system with the son's zenith angle (0* 

In the earlier investigations of radiation balance by satellites (References 14, 16, 40), the 
earth -atmosphere system was assumed to be a Lambert reflector* The albedo (References 18, 41) 
was then obtained by multiplying p* in Equation 3 by n. 


GEOGRAPHIC DISTRIBUTIONS OF RESULTS 

Geographic distributions of the outgoing long-wave radiation, albedo, and radiation balance are 
shown in Mercator projections, 70 r: N to 60°S, and polar-stereographic projections for both polar 
regions (Appendix E). 

The albedo maps for all five periods show essentially the same patterns. Albedos of more than 
55 percent were obtained over the arctic icecap and over snow -covered surfaces, especially in 
May. From May to July the mean albedo of the polar ice shield diminishes from about 75 to 60 
percent, because of the melting of the snow cover and some pack ice* The albedo over the Green- 
land ice -shield remains nearly constant at 75 to 80 percent. Although over the Arctic a mean cloud 
coverage of more than 6/10 (see Appendix F) was reported, the pattern of albedo clearly indicates 
the change of the mean pack-ice boundary with the season, demonstrating the influence of the re- 
flection properties of the earth’s surface on the albedo of the earth- atmosphere system above clouds. 
The pack ice was an average albedo of 50 to 60 percent contributes much to the flux of reflected 
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solar radiation above clouds. For comparison* the albedo over the open oceans is less than 10 
percent. (Bee Reference 21 for a more comprehensive description of radiation balance and its 
components.) 

Toward the equator the albedo decreases rapidly with decreasing cloudiness* showing minima 
of less than 20 percent over the subtropical oceans on both sides of the equator. High cloudiness* 
especially over Central Africa, and the Monsoon over Southeast and South Asia, cause albedos of 
more than 20—30 percent and 40 percent, respectively. The Sahara and Arabian deserts have al- 
bedos upward of more than 30 to 35 percent. South of 30°S, the albedo increases nearly tonally to- 
ward 60°S because of the increasing cloud cover and decreasing sun elevation. The sun did not il- 
luminate latitudes poleward of 67— 71*S during Nimbus n observations. 

The pattern in the maps of the outgoing long-wave radiation is narrowly correlated to the 
temperature of underlying surfaces (Reference 42). Very low values are obtained, especially over 
the South Pole (- winter pole), Greenland, and the ice shield over the East Siberian Sea. Else- 
where, at mean and low latitudes, fluxes of less than 0.36 cal cm“ a min' 1 are associated with cloud 
fields in areas of high cyclonic activity and areas of the intertropical convergence zone (ITC). The 
ITC is more pronounced in outgoing long-wave radiation than in the albedo, possibly because of the 
large temperature difference between high tropical cloud surfaces and the adjacent nearly cloud- 
free and warm subtropical ocean and land surfaces (these are the regions of highest loss of long- 
wave radiation to space). The temperatures over the Sahara and Arabia increase markedly from 
May to July. 

The maps of the radiation balance or the net radiation flux at the top of the atmosphere clearly 
show major regions of gain and deficit of radiative energy between 16 May and 28 July 1968. Owing 
to the suiTs declination during this season, the northern hemisphere south of 70°— 75°N and the 
southern tropics absorb more radiation than they re-emit to space. Maxima of energy gain occur 
over the very low-reflecting northern subtropical oceans, hi contrast to these, the large desert 
areas of North Africa and Arabia (both are highly reflecting and emitting areas) have only a small 
radiation surplus (in May) or even a small deficit. The latter was predicted already by Budyko 
(Reference 6) for the annual average. 

The Arctic, especially the inner ice-shield over Greenland, reflects and emits more energy to 
space than it gains by absorption of incoming solar radiation. The radiation balance over central 
arctic regions reaches positive values only during the first half of July. The southern hemisphere 
approximately south of 5°— 10°S is an area of radiation deficit during this season. Because of the 
continuous decrease of incoming solar radiation to 70°S (where the influx is almost zero) and de- 
creasing temperatures toward the Antarctic, the maximum deficit occurs in a zonal belt between 
about 60° and 70 °S. Over most land surfaces, especially the southern hemisphere, the radiation 
deficit is larger and the gam is smaller than over ocean surfaces at the same latitude. 

In all maps, the patterns of isolines over the southern hemisphere is almost zonal by latitude, 
because of the preponderance of ocean. It is very pronounced in the maps of the radiation balance, 
possibly because the meridional gradient of incoming solar radiation is very strong over the southern 
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hemisphere daring this season. Over the northern hemisphere the wide continental surfaces dis- 
solve the zones in several maxima and minima. 


GLOBAL AND HEMISPHERICAL BUDGET 

Table 1 summarizes the global and hemispherical averages of all components required in the 
balance equation; also the average albedos. The global average of the incoming flux of solar radi- 
ation s ranges between 0,484 and 0.488 cal cm’ 2 min’ 1 , changing its value according to changes in 
sun-earth distance. During the entire period the southern hemisphere receives only half as much 
solar radiation as the northern hemisphere. The global albedo was only 29,1—30.6 percent, which 
is considerably less than earlier accepted values of 33 percent and more, which were estimated or 
obtained from calculations with climatological data (References 1, 8, 43—46) and from observations 
of earth-shine on the moon (References 47, 48), Changes of the planetary value and of the 

Table 1 


Hemispherical and Global Averages of the Incoming Flux (S) and Reflected (R) Flux of Solar Radiation 
of the Albedo (R/S), of the Outgoing Long-Wave Radiation (E) and its Equivalent Black- Body Temperature, 
of the Absorbed Solar Radiation (S^R), and of the Net Radiation Flux (Q). 




S 

R 


E 


S-R 

Q 

Period 

(1966) 

Hemisphere(s) 

cal cm" 

2 min” 1 

K/ a 
x 100 
(percent) 

Equivalent 

Temperature 

m 

cal cm“ 2 min" 1 


Northern 

0.659 

0,216 

32.8 

255.2 

0.345 

0.443 

+ 0,098 

16-31 May 

Southern 

0.316 

0.078 

24.7 

252.7 

0.333 

0.238 

-0.095 


Globe 

0.488 

0.147 

30.1 

254.0 

0.339 

0.341 

+0.002 


Northern 

0.674 

0.221 

32.8 

256.0 

0.350 

0.453 

+0.103 

1-15 June 

Southern 

0.297 

0.076 

25.6 

253.2 

0.334 

0,221 

-0.113 


Globe 

0.486 

0.149 

30.6 

254.8 

0.342 

0,337 

-0,005 


Northern 

0.676 

0.217 

32.1 

256.5 

0,353 

0,459 

+0,106 

16-30 June 

Southern 

0.291 

0.074 

25.4 

253.5 

0.336 

0.217 i 

-0.119 


Globe 

0.484 

0.146 

30.1 

255,2 

0.345 

0.338 

-0.007 


Northern 

0.670 

0.210 

31.4 

257,0 

0.356 

0,460 

+ 0.104 

1-15 July 

Southern 

0.298 

0.072 

24.1 

253.5 

0.336 

0.226 

-0.100 


Globe 

0,484 

0.141 

29.1 

255,2 

0.346 

0.343 

-0,003 


Northern 

0.653 

0.203 

31.1 

256,8 

0.355 

0,451 

+0,096 

16-28 July 

Southern 

0,317 

0.084 

26.5 

253,2 

0.334 

0.233 

-0,101 


Globe 

0.485 

0.143 

29,5 

255,2 

0.345 

0,342 

-0.003 
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hemispheric averages might be caused more by insufficient data coverage over parts of the globe 
than by albedo changes. The data coverage of the globe throughout all five subperiods was not al- 
ways complete, as indicated by stippled fracture and dashed iso lines in the maps. Thus, all av- 
erages in Table 1 might be slightly biased by lack of data in these gaps. The average albedo is 
rated higher over the northern than over the southern hemisphere, since during this season there 
were no measurements obtained over die Antarctic, Further, the southern hemisphere is predom- 
inantly covered with low- reflecting ocean surfaces. 


From values of the outgoing flux, the mean "effective” temperature of the planet earth was 
found to be between 254 and 255.2°K, which is considerably higher than the 251 °K previously ac- 
cepted assuming radiative equilibrium for the earth- atmosphere system and a 35-percent global 
albedo. Here, the northern hemisphere is about 3 degrees warmer than the southern. From May 
to July, the emitted long- wave radiation rises slightly in the northern and decreases slightly in the 
southern hemisphere. 


Figure 5 shows separately the temporal 
trends of the semihemispherical and global 
averages of the albedo and of the outgoing 
long- wave radiation. Over the northern hem- 
isphere the outgoing flux rises slightly while 
the albedo simultaneously decreases. Both 
trends are caused by the warming of the north- 
ern hemisphere during this season and the 
melting of parts of the arctic ice cap. There 
seems to foe no similar trend in the opposite 
direction over the southern hemisphere. The 
albedo, there, was found to be extremely low 
(24.1 percent) for the first half of July; but no 
explanation has been found. The period of 
available Nimbus II measurements was too 
short to show definite seasonal trends in these 
large-scale averages of the albedo and of the 
outgoing long-wave radiation. 
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Figure 5 — Hemispherical and global averages of the 
outgoing long^vave radiation, flux (E) and the albedo 
(R/S), 


The last column of Table 1 lists the re- 
su its on glob al r adi ati on b alanc e , Th ey i ndic ate 

for June and July 1966 a minimal radiation deficit of the planet earth, which is less than 2 percent 
of both components— the absorbed solar radiation and the emitted long- wave radiation. For the 
second half of May a slight radiation gain was found. Since the value of the solar constant is known 
at present only within ±2 percent, no final conclusions can be made on the real budget of the earth- 
atmosphere system during available Nimbus H measurements. However, a similar slight deficit 
found for this season by Simpson (Reference 3) suggests that the earth may have a slight radiation 
deficit during this season, possibly due to the larger earth- sun distance. Further measurements 
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over periods of 1 or 2 years (including simultaneous measurements of the solar constant) will lead 
to more conclusive results on the earth radiation budget. 

The maps in Appendix E also show that the northern hemisphere permanently gains in radia- 
tive energy , while the southern hemisphere has a radiation deficit But the effective emission tem- 
perature (T) for both hemispheres at the top of the atmosphere indicate that one or the other is not 
heating or cooling definitely (Table 1), Thus, atmospheric circulation as well as energy transport 
by ocean currents accounts for the energy transport from regions of radiation surplus to those of 
deficit. This, possibly, should also act southward through the equator, although it does not appear 
pronounced in mean streamline patterns for this season (Reference 49). 


Spatial and temporal changes of the radiation budget over different latitudinal zones are better 
shown in the zonal averages plotted in Figure 6: incoming solar radiation, albedo, outgoing long- 
wave radiation, and radiation balance. The incoming flux of solar radiation changes considerably 
north of 40°N and over the southern hemisphere. These changes are explained by changes of the 
sun*s declination which is greatest at the summer solstice. 

The zonal averages for the outgoing long-wave radiation show the trend expected over both 
hemispheres. From May to July the northern hemisphere becomes slightly warmer north of 30^N, 
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Figure 6— Zonal averages of (a) the incoming solar radiation, (b) the outgoing long -wave radiation, 
(c) the albedo, and (d) the radiation balance; at the top of the atmosphere. 


but slightly colder south of 30 because of northward shift of subtropical anticyclones and of the 
ITC (and all cloud fields associated with it), , The thermal gradient from the equator to the North 
Pole decreases slightly, which may result in a sii^it decrease of the circulation (Figure 27 of Ref- 
erence 8), The opposite behavior is observed over the southern hemisphere, where warming in the 
southern tropics (caused by the northward shift of the ITC) and cooling in the Antarctic cause the 
thermal gradient (and possibly the circulation) to increase slightly. Toward the end of July, radia- 
tive cooling at the outer belt of the Antarctic cannot be compensated for completely by advection. 

The slight heating over the northern hemisphere is accompanied by a decrease of the albedo 
north of 60°N, which is due to melting of the snow cover over the arctic pack-ice, and of parts of 
itself. The radiation balance, however, was found positive over the arctic only during the first 
half of July, But the deficit is small enough, for parts of the energy transported to the arctic by 
advection and ocean currents to be available to melt parts of the ice and increase the mean 
temperature there. 

From 50°N to 20°S the albedo shows no definite temporal changes associated with the north- 
ward movement of the ITC, Remarkably lower albedos occur during the first half of July over the 
southern hemisphere; a fact not yet explained. The radiation balance of the earth- atmosphere sys- 
tem south of the equator to 50 changes its value almost entirely because of solar -influx changes. 
The maximum deficit of more than 0,24 cal cm' 2 min" 1 is observed at the summer solstice at 60 °B 
and decreases there rapidly toward the end of July because of decreasing temperature and increas- 
ing insolation, A more detailed description of the results obtained over both polar caps is given 
elsewhere {Reference 21), 

For Nimbus H, even the semimonthly averages might not very well represent the mean radia- 
tion balance of the earth- atmosphere system over areas with pronounced diurnal cycles of cloud- 
iness and temperature and small transient disturbances; such was the character of the orbit. These 
areas are located mainly in tropical and subtropical regions. 


COMPARISON OF THE RESULTS WITH DIRECT MEASUREMENTS 

A comparison of our averaged results with direct (and ideally simultaneous) measurements 
should clearly reveal errors inherent in our methods. Because of a satellite’s height such meas- 
urements should be made from high-flying airplanes or from balloons over areas whose radiation 
fields are nearly homogeneous in time and in space. 

Such areas occur over those regions of the earth where synoptic disturbances (e,g, changes in 
cloud cover and cloud height) cause— in shorter time intervals— only small variations in both fields, 
in the outgoing long- wave radiation and in reflected solar radiation (or albedo). They have been 
found by computations of standard deviations of the albedo and of the outgoing long- wave radiation 
of all single daily averages within a grid field. Figure 7 shows results for the albedo during June, 
1966, The standard deviation was less than 5 percent of the average value for both quantities, the 
albedo and the outgoing radiation flux {not shown here) over the Arctic (especially Greenland) and 
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over the Sahara* indicating the predominant influence of surface conditions on the outgoing radiation 
in these areas* Standard deviation of the albedo over the subtropical oceans is less than 5 percent; 
but actually this is large, considering that the total mean albedo there is 15 percent or less. Over 
land areas at the same latitude the standard deviation is small and the mean albedo is much higher. 
Also, the standard deviation of the outgoing long-wave radiation is very low over the Antarctic, 

In- situ measurements of outgoing long-wave radiation were recently discussed by Kuhn*, who 
conducted radiometer- sonde measurements over six antarctic stations during the same period {16 
May— 28 July 1966), His results agreed with our results within ±2 percent. It is now highly desirable 
to make a similar comparison over much warmer areas, such as the Sahara or Australian deserts. 

No simultaneous measurements of the outgoing shortwave radiation were available for com- 
parisons with the albedo values obtained from the Nimbus H measurements. But only a few authors 
reported on results which they obtained on the directional reflectance of surfaces during the same 
season although in different years. Kasten {Reference 50) obtained over Greenland values of 
around 72 percent under clear skies in June; this compares quite well with our albedos between 70 
and 80 percent which are affected by the atmosphere. The cloud cover over Greenland was about 
4/10 to 6/10 during June (see Figures F2 and F3), Diamond and Gerdei (Reference 51) reported 
similar values; they obtained 77 percent in the forenoon and 87 percent in the afternoon under a 


* Kuhn, P. H . , "Nimbus II and Balloon Radiation Analyse s, " to be published in /. Atm. Sri. , private communication, 1967, 
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clear sky, Hubley (Reference 52) obtained over the Lamon Creek Glacier in Alaska a pronounced 
variation of the directional reflectance with local time. His values range between 80 and 65 per- 
cent at forenoon and 65 to 85 percent in the afternoon, Predoehl and Spano (Reference 53) and 
Spano (Reference 54) reported on airplane measurements over the Antarctic ice cap, "These meas- 
urements, though conducted during a different season, show by contrast the extent of our results 
over the Arctic regions. Over the old winter ice they found values between 65 and 14 percent; new 
ice had considerably less directional reflectance. These values seem to confirm the high albedos 
of 70 to 80 percent, which were obtained from the Nimbus n measurements over Greenland by the 
computational methods discussed above. The straight assumption of a Lambert- surface led to con- 
siderably smaller values between 55 and 60 percent (Reference 41). 

Roach (Reference 55) obtained over the Sahara between Ei-Adem and Karthoum directional 
reflectances of about 32 percent. Our values over this region range between 30 and 40 percent and 
may be affected largely by the cloudiness. This rough comparison shows that our computational 
results obtained from Nimbus H radiance measurements agree quite well with other measurements. 
Definite conclusions on the radiation balance of the planet earth, and for specific areas, require a 
more reliable confirmation. 

COMPARISON WITH OTHER RADIATION BALANCE STUDIES 

From climatological data Baur and Philipps (Reference 4) and London (Reference 1) obtained 
results for die same season (summer) for the northern hemisphere only, Gabites (Reference 10) 
computed monthly averages of the radiation balance and its components only south of 40 °S, using 
observations of the International Geophysical Year 1957 — 1958. Only Simpson (Reference 3) and 
Vinnikov (Reference 7) obtained results on a global scale. Houghton (Reference 5) and Budyko {Ref- 
erence 6) reported only on annual averages. 

Measurements of previous TIROS satellites were evaluated by several authors: Rasool and 
Prabhakara (Reference 16), Winston (Reference 17), Bandeen, et al, (Reference 14). The outgoing 
flux of long- wave radiation was computed from measurements of infrared radiation between 8 and 
12 microns, since measurements in a wide spectral range were unavailable (on TIROS IV) or very 
inaccurate (on TIROS V3I), These measurements and those of the outgoing solar radiation were 
highly affected by instrumental degradation. The results, then, may still contain some errors, al- 
though careful correction methods were applied (Reference 14), The albedo of the earth-atmosphere 
system was computed in these earlier studies assuming complete isotropy of the reflection proper- 
ties of the earth- atmosphere system. Although the original measurements were corrected for deg- 
radation, the albedos were still too low. Arking and Levine (Reference 56) applied an angular cor- 
rection on TIROS VEt measurements. But they had to use a factor of 1.56 to obtain a 32-percent 
global albedo; this albedo value was necessary to balance in the annual average the incoming solar 
radiation with the outgoing radiation fluxes (Reference 14), From Explorer VH data, House (Refer- 
ence 15) obtained averages for December 1959 through May 1960 only. Thus, our results cannot be 
compared with those from earlier satellite measurements. 

Figure 8 summarizes zonal averages of earlier investigations and compares them with those 
obtained from Nimbus II measurements throughout July, The estimates of the albedo and outgoing 
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Figure 8— Comparison of zonal averages of (a) the albedo, (b) the outgoing long -wave radiation, and (c) the 
radiation balance, obtained from Nimbus IE measurements during the period 1—28 July 1966, with results of 
other authors* 
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long -wave radiation from climatological data resulted in higher albedo values and lower outgoing 
radiation values between 50°N and 40 °S than were obtained from Nimbus H data. This may be due 
to overestimating the cloud cover, the cloud height, and die reflection properties of clouds. A com- 
parison of die albedos with the surface albedo, as shown in part (a) of the figure, may give some 
idea of the contribution of clouds to the albedo. (Between TQTST and 90^ Posey and Clapp (Refer- 
ence 57) may have considerably overestimated the albedo of the polar ice cap and surrounding 
areas, not taking into account the open sea, which has considerably less albedo,)* 

All authors found the radiation gain and deficit lower between 50°N and 50°S-see Figure 8(c). At 
higher northern latitudes only London *s (Reference 1) results show a good agreement. Over the cen- 
tral Arctic (where Nimbus H measurements led to a slight radiation gain of 0.001 cal cm" 2 min" 1 for 
the monthly average), Vowinckel and Orvig (Reference 9) found a much higher surplus of +0.05 cal 
cm -2 min" 1 , while Fletcher (Reference 8) obtained a deficit of -0,02 cal cm" 2 min“ l , Gabites T 
results (Reference 10) over high southern latitudes show the closest agreement with our results, 
Simpson (Reference 3) had estimated too high surface temperatures over the Antarctic continent. 


CONCLUSIONS 

The main object of these investigations has been to calculate— from Nimbus 13 measurements 
of the radiance of outgoing long- wave and short-wave radiation— the outgoing flux of long-wave ra- 
diation and the albedo of the earth- atmosphere system. Both of these quantifies and the incoming 
flux of solar radiation (obtained on the basis of the solar constant S 0 = 2,0 cal cm” 2 min" 1 ) were 
used to compute the radiation balance, which is essentially the net radiation flux at the top of the 
atmosphere. Results are presented in geographic maps for the entire globe for five semimonthly 
periods covering the entire period of available Nimbus II measurements (16 May to 28 July, I960), 
The spatial resolution was chosen to be between 500 km X 500 km (at the equator) and about 280 km 
x 280 km {poleward of latitude 60°), These maps show for the first time global distributions of the 
albedo, the outgoing long-wave radiation and the radiation balance as obtained from satellite 
measurements. 

Other main results are: 

1. Consideration of the anisotropic reflection characteristics of the earth- atmosphere system, 
in simple models, in albedo computations led to results that reasonably agree with meas- 
urements by other authors. But only a few comparisons were possible, 

2, Outgoing long-wave radiation-flux measurements over the Antarctic agreed within 2 per- 
cent with simultaneous radiometersonde measurements. Similar comparisons over warmer 
areas are desired. 


also the orbit [tl characteristics of Nimbus II have caused errors over areas where the diurnal cycle of cloudiness and cloud height 
(and thus of the outgoing long- wave radiation and the albedo) is very pronounced as compared with day-to-day changes due ro transient 
disturbances* Merritt, et al, (Reference 58) , for instance, obtained a semi-diurnal cycle of cloudiness over Indonesia and Malaysia, with 
maxima in the early morning and late afternoon, while at local noon the cloudiness was somewhat less. Nimbus II daylight measure- 
ments ate obtained over these areas around local noon; thus, the albedo may have been somewhat lower than if measured throughout the 
whole day* But the differences between the older studies and our investigations are too pronounced to be due only to these diurnal 
effects. 
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3. Global averages of the radiation balance of the earth-atmosphere system range, for all 
five subperiods, within 0.002 cal crrT 2 min' 1 gain and 0*007 cal cm" 2 miiT 1 deficit {which 
is less than 2 percent of both components, the absorbed solar radiation remaining at the 
planet earth, and the emitted long- wave radiation being lost to space). These results indi- 
cate that the earth-atmosphere system is— even in smaller time periods— close to radiative 
equilibrium, {But, the results might be biased also by a systematic error in the solar 
constant*) 

4, Zonal and global averages of the albedo and outgoing long- wave radiation are lower and 
higher, respectively, than those obtained from climatological data. This disagreement 
occurs mainly at low latitudes. This suggests that ground-based observations considerably 
overestimate the cloud cover. 

For future efforts (recommended in research programs, such as the "Global Atmospheric Re- 
search Program," Reference 62 ? and the "World Weather Watch," Reference 63) the directions are 
very clear: 

1. The solar constant must be measured simultaneously with other measurements of radiation 
emerging from the earth. 

2. Higher spatial and temporal resolutions of the results require more sophisticated computa- 
tional efforts, which permit the study of climatological peculiarities over different areas, 

3. Simultaneous but independent measurements of the outgoing long- wave radiation flux and 
reflected solar radiation must be conducted over selected areas, to check the quality of the 
results obtained from satellite measurements of the radiance. 

The results obtained so far from Nimbus II measurements on the radiation balance of the earth- 
atmosphere system clearly demonstrated the validity of satellite measurements. Results over 
several years {as well as studies of their relation to the general circulation pattern) have a funda- 
mental bearing on the energy exchange of the planet earth. Once approved, similar methods could 
be applied for other planets. 
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Appendix A 

Computation of Angles e, : 


Solar Zenith Angle: 


£ 


COS 


j^cos 4 > cos S cos - X ) + s in $ s in si 


(Al) 


where X, 4 > are the geographic longitude and latitude, respectively, of an observed surface element 
(\ is counted westward; 4> is positive north and negative south of the equator), X Q is the Green- 
wich hour angle, and s is the declination of the sun. 


Zenith Angle of Measurement 8: 


6 


sin \ — ^ — sin a) 

sin^ 1 (K * sin a) , 


(A2) 


where a is the nadir angle of measurement, which can be obtained from the data tape (Reference 22), 
The earth's mean radius R is assumed to be 6311 km and the satellite height H to be constant; 

H ~ 1140 km. 


Then K = (R+h)/R = 1,17893, 


Azimuth Angle of Measurement \p: 

The azimuth angle $ is deter mined with respect to the ray of incident solar radiation (Fig- 
ure 1), by 


_ , / cos r- cos £ cos 8 \ 
0 - cos \ sin i sin 6 } 


(A3) 


where 


cos r - 


- cos cos (X g — Xq ) cos <£ s cos (x - X Q ) j cos £ + [k sin 4> s * sin <£>] sin £ 
|/k 3 + 1 - 2K |cos <p s cos ^cos (X - X s ) + sin$ s sin$j 


(A4) 


where X g and <£ s are the geographic longitude and latitude, respectively, of the subsatellite point, 

26 


The Length of a Day AX : 

Sunrise and sunset are defined as the instants when the upper edge of the disk of the sun is on 
the horizon at normal refraction {Table IX in Reference 39), Assuming that a semidiameter of the 
sun subtends a 16 “minute angle , with a constant 34“ minute refraction (Reference 59), then the 
length of a day AX* {in degrees of arc on the celestial sphere) is 


AX* = 2 


AX 


-Q. 0143 - sin sin S 
cos 4 > cos 0 


2 cos" 1 ( arg) , 


(A5) 


where A\ is 180° (= 12 hours), if (arg) 5 - 1. 



Appendix B 


Computation of Outgoing Long-Wave Radiation Flux E 

II N f (#) is the measured filtered radiance in the spectral range from 5.0 to 30.0 microns, then 
the unfiltered radiance N {Figure 2) is:* 

N(tf) = ^ [ok 0439 + 5.4318 x 10" 3 ■ N f (£) + 1.0245 *10~ s * N f (0) 2 

+ 0 3 * { 1. 36 x 10“ 8 - 7.9262 x 10“ 10 ■ N f (0) + 1 . 1483 x 10^ 11 ■ N f (0) 2 )J , (gi) 


where N(£) and N f ( 6 ) are in cal cm' 1 min 1 sr ~ 

The outgoing flux E, then* is determined by 

E " 2.901 - W{6)/[l + 1.933 x!0^ 4 ■ 0 -4.247 x 10" s ■ 6 2 

+ 6. 149 x 10' 7 ■ £ 3 - 7.807 x 10’ 9 ■ & 4 ] cal cm -3 min' 1 . (B2) 

The denominator in Equation B2 was found statistically by Lienesch* from measurements of TIROS 
satellites. It describes the zenith- angle dependence of outgoing long-wave radiation. 


•t.ienesch, J. H., private communication, 1966. 
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Appendix C 

Relations to Determine Directional Reflectance r 
from Bidirectional Reflectance p 


General Equation for HO: 

To integrate Equation 5 we need to know how the directional reflectance of an observed surface 
element varies with the solar zenith angle. Since only one or two satellite measurements were pos- 
sible over the same area throughout the day, an empirical model was set up, to determine r(£) from 
the bidirectional reflectance p' (8, £ ' } , as obtained from a measured radiance (Equation 3): 


r(0 


*( 0 _ 
r(£ = 0) 


r(t = _1> . '(O 

r(t') P{8, 4>, £') 


”p' (8, i p, V) , 


(Cl) 


(See Figures 3 and 4) 

where £' and £ are the solar zenith angle at the moment of measurement and at any arbitrary time, 
respectively, throughout the day. 


Derivation of the Relation r(£)/r (£ = 0): 

The general relation for the change of the directional reflectance with the sun's zenith angle 
was derived from several measurements. A main condition for the selection of measured reflec- 
tances was: that these measurements should be made in a wide spectral interval comparable with 
that of the Nimbus n MRIR (0.2— 4.0 microns). Further, these measurements should have been 
obtained at a certain distance above reflecting surfaces (ground or clouds), so that they should 
also contain the effect of the atmosphere above those surfaces. 

Up to that time (when these calculations were made), only a few measurements were available; 
these were obtained from airplanes (References 31, 32) and from balloons (Reference 30). The 
underlying surfaces were clear ocean, clouds, and snow. These measurements were obtained only 
in three narrow ranges of the solar zenith angle; 10° £.£ £■ 20°; 50° £ £ £ 60"; 65° £ £ £-75°. Arking (Ref- 
erence 33) obtained results on the angular behavior of the reflectance from a statistical analysis 
of TIROS IV measurements. 

The heterogeneity of these available data required the assumption of a multiplicative law, to 
give an equation for the mean dependence of r on £: 


r(£) _ Pi 0, 0, £) . r(£) _ p(0, 0, 0) , 

r(£ = 0) P<0, 0. 0) p(0, 0. £) r(£ = 0) 

where ^(0, o, £) = v • p' (0, 0, £) at the subobserver point on the earth's surface. 


(C 2) 
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Figure Cl shows the dependence of the product of the last two factors in Equation C2 on £ , 
which describes the change of the ratio between the directional reflectance and the bidirec- 
tional reflectance at the subobserver point with the sun T s zenith angle, A smooth curve drawn 
through the various kinds of data point was adopted in these computations. For comparison, the 
results of Arising* and Levine f are drawn in Figure Cl. 

The factor p(Q f o, O/p(0, 0, 0) , which describes the dependence of the bidirectional reflectance 
of the subobserver point on the solar zenith angle, was obtained from Arking’s and Levine's statis- 
tical analysis of TIROS W measurements. Figure C2(a) shows those authors' results. The rela- 
tion of p(0, 0, l)fp{ 0, 0,0) to l adopted for this work was obtained by curve -smoothing and is shown 
in Figure C2(b). This curve has some similarity with scattering functions obtained theoretically 
for water droplets (see for instance, Reference 61); this may be incidental and be due to data 
sampling. 


The final relation of r (£}/r{£ = 0) to £ was obtained through Equation C2 and is shown in Fig- 
ure 4 in the main body of this report. 


Diagrams X = r/p: 


According to the third factor in Equation Cl, we need relations to convert the bidirectional 
reflectance p 1 obtained from Equation 3 to the directional reflectance r. These relations, assumed 


to represent the mean behavior of the earth- atmosphere system, replace the measurements which, 
in fact, are required under all angles 9, \p t £ over one area. Bartman (Reference 30) and SaXomonson 
(Reference 32) obtained measurements over different surfaces for many values of 9 and 0. Cherrix 


CJ 

II 


K0 

0 


-X LEVINE, (TIROS Vlh 0.55^ -0.75 P ) r 
PRIVATE COMMUNICATION, 1967 
— ARKING, (TIROS IV z 0.2^ -6.0^ ), 
PRIVATE COMMUNICATION, |96? 

A REFERENCE 30 ( CLOUDS ) 
o REFERENCE 29 (CLOUDS) 

© REFERENCE 31 (CLOUDS) 
o REFERENCE 30 ( CLEAR OCEAN) 
w REFERENCE 29 (SNOW) 

REFERENCE 60 (OCEAN) 

• REFERENCE 32 



10 20 30 40 50 60 70 80 

SOLAR ZENITH ANGLE <( degrees ) 


and Sparkman (Reference 31 at the time had re- 
sults for azimuthal angles ^of0°, 180°, 45°, 135°, 
90° only; thus, in order to provide models at least 


G^-O TIROS VII (0.55u-0-75p) x-x ABOVE MEAN 



SOLAR ZENITH ANGLE f (degrees) 


Figured— The relative change of the ratio between the 
directional reflectance (r) and the bidirectional reflec- 
tance^ = pfn) at the sybobserver point {9 - 0° f \p ~ 0°) 
with the solar zenith angle {£). 


Figure C2— The relative change of the bidirectional re- 
flectance ( p ' = p/tt) at the subobserver point (9 ~ 0° f 
0 = 0°) with the solar zenith angle {£); experimental 
and smoothed curves. 


* Arlclngj A., private communication, I?67* 

'Levine, J. S., "The Planetary Albedo Based on Satellite Measurements Taking into Account the Anisotropic Nature of the Reflected 
and Backscattered Solar Radiation,” Master's Thesis, New York Univ., Graduate School of Arts arid Science (private communication). 
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on the average behavior, all available data have been summarized only lor these azimuthal angles* 
Symmetry was assumed with respect to the principle plane (^: 0°— 180°), 


Figures 03 through 05 show the mean curves r/p , which were also obtained by smoothing. 
The single measurements, average lor - 10 degrees, are indicated by several signatures ac- 


o STRATOCUMULUS (REFERENCE 31) 

O STRATOCUMULUS (REFERENCE 30 FLIGHT 35) 

• STRATOCUMULUS ( REFERENCE 30 FLIGHT 20). 



ZENITH ANGLE OF MEASUREMENT (degrees) 

Figure C3— The ratio between the directional reflec- 
tance (r) and the bidirectional reflectance^' = pfn )af 
values of azimuthal angle {pi 0®, 180°, 45°, 135°, 90°) 
as obtained From measurements at solar zenith angles 
( 10 ®<£< 20 °)* 


□ SNOW (REFERENCE 29) 
o STRATOCUMULUS (REFERENCE 32} 



ZENITH ANGLE OF MEASUREMENT ( degrees ) 


Figure C4— The ratio between Hie directional reflec- 
tance (r) and the bidirectional reflectance (p* “ p/n) at 
values of azimuthal angle {p: 0°, 180°, 45®, 135°, 90°) 
as obtained from measurements at solar zenith angles 

(50® <1 160°) * 


cording to the data source* AIL values r/p are 
normalized with respect to the relation shown in 
Figure Cl* To draw the mean curves in Fig- 
ures C3 through C5, most weight was given to 
measurements which were obtained in all five 
planes = 0°, 180°, 45°, 135°, 90°) over cloud 
surfaces. 
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O STRATOCUMULUS (REFERENCE 31 ) 
o SNOW (REFERENCE 29) 

O STRATOCUMULUS (REFERENCE 30, FLIGHT 36) 
• STRATOCUMULUS ( REFERENCE 30, FLIGHT 37) 




ZENITH ANGLE OF MEASUREMENT (degrees) 


Figure C5*~The ratio between the directional reflec- 
tance (r) and the bid? reef I anal reflectance ( p * = pfi r)qt 
values of azimuthal angle {pi 0°, 180°, 45°, 135°, W 5 ) 
as obtained from measurements at solar zenith angles 

(65° < £< 75°). 
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The deviations of all single values from the mean smoothed curves in Figures C3— C5 is 
rather large, showing that the reflection properties of single cloud decks and snow covers deviate 
from each other- But /\e smooth curves indicate, at least, a general trend, which is used here to 
account for the angular dependence of Nimbus H measurements of reflected solar radiation. 

The diagrams shown in Figures 3, C6, and C7 were obtained from the mean curves in Fig- 
ures 03 through C5 by graphical interpolation. The angular dependence of the reflectance observed 
for 10° < £ < 20 ° was adopted for all Nimbus H measurements for 0° < £ 5 35"; that observed for 
50° < £ < 60 for 35 c < £ 1 60°; and that observed at 65° < £ < 75°, for 60° * £ 1 80°. The sparsity of 
observations then available made it impossible to determine more detailed models. 


0 = 180 ° 



Figure C6 —Dependence of the ratio X - v/p on the 
angles# and 0 of the measurement at very high sun 

{ 0 ° Si <35% 


&= I 80 a 



Figure C7 — Dependence of the ratio X - v/p on the 
angles 8 and \p of the measurement at solar zenith angles 
( 35 ° < £ 1 60 °) . 
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Appendix D 


Summary of Computer Program 


Figure D1 is a system flow chart that summarizes all phases in the computational procedure. 

Phase 1 computes all required radiation fluxes for all orbits of one day or parts of a day. The 
NMRT— MRIR tape (see Reference 23) contains the measurements. Control parameters are re- 
quired to read the tapes as well as those parameters required in several equations. The results 
of each orbit for one day or parts of a day will be stored in matrices (which correspond to stand- 
ardized matrices) for geographic mapping of the satellite data. 


FOR TAPE, THESE DATA SETS 
MAY BE WRITTEN AS SEPA- 
RATE FILES ON THE SAME 
TAPE. THEY ARE WRITTEN IN 
THE ORDER: N. ROLE, MERCATOR, 
S. POLE. 


NOTE: ANY SEQUENTIAL DATA 
STORAGE DEVICE CAN BE 
USED INSTEAD OF THE TAPES 
SHOWN ON THIS CHART. 



PHASE 2 & 3 LOGIC SIMILAR 
TO THAT SHOWN FOR NORTH POLE. 


THESE SEPARATE DATA SETS 
ARE ACTUALLY READ BY 
PHASE 3 AS ONE DATA SET. 
THEY MAY BE WRITTEN AS 
ONE DATA SET (USING THE 
DISP = MOD OPTION OF JCL) 
OR THEY MAY BE CONCAT - 
ENATED ON INPUT. 


Figure Dl— System flow chart of the computer program. 
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Phase 2 summarizes all matrices and gives averages within each gridfield for each day (24- 
hour interval, counted from 00,00 to 24,00 GMT) separately. 

Phase 3 computes averages for a number of days, as given by the control parameters. Zonal 
averages will be determined only for the Mercator matrices between 75^ and 75°S, Separate pro- 
grams were set up to compute standard deviations, relative dispersions, and differences between 
day and night values of the outgoing long-wave radiation for each gridfield. Furthermore, zonal 
averages can be determined for several zones on the globe. Maximum and minimum values of the 
albedo and the outgoing long -wave radiation (day and night separately) are determined for several 
climatological studies. 
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Appendix E 

Maps of Results 


The geographic distributions of the albedo of the outgoing long-wave radiation and of the ra- 
diation balance or net radiation flux are shown here, separately, in maps for all five semimonthly 
periods* (See Figures El through E40.) 

In the polar sfcerographic projections, the spatial resolution for each gridfield ranges from about 
280 km X 280 km at the pole to about 260 km x 280 km at 60° latitude* In the Mercator projections, 
a gridfield at the equator covers an area of 5* X 5° of longitude and latitude {500 km x 600 km) and 
at 60° latitude an area of 5° x 23° of longitude and latitude (280 km x 280 km)* Isolmes were drawn 
by hand almost everywhere at equidistant intervals* Intermediate isolmes were added over special 
areas (e.g*, Sahara, subtropical oceans, and both polar caps) to demonstrate a specific pattern* 
Dashed isolines were obtained by interpolation and extrapolation over areas where measurements 
were made on only one or two days (or not at all) during a semimonthly period. Larger gaps in the 
data coverage are marked with stippled fracture. Such gaps were found for daytime measurements 
between about 110°E and 170°E of longitude and for night measurements between about 30°W and 
10°E* An orientation error was found in the Mercator projections after the maps were drawn; all 
data were misplaced by 2*5 degrees of longitude eastward. 
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Figure E6— Albedo of the earth -atmosphere system between 70°N and 60°$ during the period 16 — 31 May 1966, 









Figure E 10— Albedo of the earth -atmosphere 
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between 70 & N and 60°S during the period 16—28 July 1966, 
















Figure 15— Outgoing long- wave radiation flux of the top of the atmosphere 
over the northern hemisphere during the period 16—28 July 1966* 











Figure El 8— Outgoing long-wave radiation flux at the fop of the atmosphere 
between 70°N and 6Q D S during the period 16—30 June 1966* 
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Figure E 19— Outgoing long^wove radiation flux at the top of the atmospbi 
between 70°N and 60°S during the period 1 —15 July 1966, 
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Figure 20-Outgoing long-wave radiation flux at the top of the atmosphere 
between 70°N and 60 e S during the period 16—28 July 1966. 







































Figure E33— Radiation balance of the earth atmosphere system between 
70°N and 6D°S during the period 16—30 June ! 966+ 




Figure E34— Radiation balance of the earth -atmosphere system between 
70°N and 60*S during the period 1—15 July 1966, 
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Figure E35— Radiation balance of the earth -atmosphere system between 
70°N and 60°S during the period 16—28 July 1966* 
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Appendix F 

Mean Cloud Cover Over the Northern Hemisphere 

Maps of the mean cloud cover (Figures FI through F5) over the northern hemisphere are 
added to this report. They were obtained by the USAF Environmental Technical Applications 
Center, Washington, D. C.,* from ground observations. 



Figure FI— Meon cloud cover over the northern hemisphere during the 
period 16— 3 T Moy 1966 from operotiond data. 


•USAF Environmental Technical Applications Center, Washington, D, C., (private communication through Maj. E. Kreins), 1967. 
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